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Abstract 



In the past few years, small scale anisotropy has become a primary focus in the 
search for source of Ultra-High Energy Cosmic Rays (UHECRs). The Akeno Giant 
Air Shower Array (AGASA) has reported the presence of clusters of event arrival di- 
rections in their highest energy data set. The High Resolution Fly's Eye (HiRes) has 
accumulated an exposure in one of its monocular eyes at energies above 10 19 5 eV 
comparable to that of AGASA. However, monocular events observed with an air 
fluorescence detector are characterized by highly asymmetric angular resolution. A 
method is developed for measuring autocorrelation with asymmetric angular resolu- 
tion. It is concluded that HiRes-I observations are consistent with no autocorrelation 
and that the sensitivity to clustering of the HiRes-I detector is comparable to that 
of the reported AGASA data set. Furthermore, we state with a 90% confidence level 
that no more than 13% of the observed HiRes-I events above 10 19 ' 5 eV could be 
sharing common arrival directions. However, because a measure of autocorrelation 
makes no assumption of the underlying astrophysical mechanism that results in 
clustering phenomena, we cannot claim that the HiRes monocular analysis and the 
AGASA analysis are inconsistent beyond a specified confidence level. 

Key words: cosmic rays, anisotropy, clustering, autocorrelation, HiRes, AGASA 
PACS: 98.70.Sa, 95.55.Vj, 96.40.Pq, 13.85.Tp 



1 Introduction 



Over the past decade, the search for sources of Ultra-High Energy Cosmic 
Rays (UHECRs) has begun to focus upon small scale anisotropy in event ar- 
rival directions. This refers to statistically significant excesses occurring at the 
scale of < 2.5°. The interest in this sort of anisotropy has largely been fueled 
by the observations of the Akeno Giant Air Shower Array (AGASA). In 1999 
[1] and again in 2001 [2], the AGASA collaboration reported observing what 
eventually became seven clusters (six "doublets" and one "triplet" ) with esti- 
mated energies above ~ 3.8 x 10 19 eV. Several attempts that have been made 
to ascertain the significance of these clusters returned chance probabilities of 
4 x 10- 6 [3] to 0.08 [4]. 

By contrast, the monocular (and stereo) analyses that have been presented by 
the High Resolution Fly's Eye (HiRes) demonstrate that the level of autocor- 
relation observed in our sample is completely consistent with that expected 
from background coincidences [5,6,7]. Any analysis of HiRes monocular data 

* Corresponding author. E-mail address: stokes@cosmic.utah.edu (B.T. Stokes) 



2 



needs to take into account that the angular resolution in monocular mode is 
highly asymmetric. 

It is very difficult to compare the results of the HiRes monocular and AGASA 
analyses. They are very different in the way that they measure autocorrelation. 
Differences in the published energy spectra of the two experiments suggest 
an energy scale difference of 30% [8,9]. Additionally, the two experiments 
observe UHECRs in very different ways. The HiRes experiment has an energy- 
dependent aperture and an exposure with a seasonal variability [8]. These 
differences make it very difficult get an intuitive grasp of what HiRes should 
see if the AGASA claim of autocorrelation is justified. In order to develop this 
sort of intuition, we apply the same analysis to both AGASA and HiRes data. 



2 The HiRes-I Monocular Data 

The data set that we consider consists of events that were included in the 
HiRes-I monocular spectrum measurement [8,10]. This set contains 52 events 
observed between May 1997 and February 2003 with measured energies greater 
than f 19 5 eV. The data set represents a cumulative exposure of ~ 3000 km 2 -sr-yr 
at 5 x I0 19 eV. This data was subject to a number of quality cuts that are 
detailed in the above-mentioned papers [8,10]. We previously verified that this 
data set was consistent with Monte Carlo predictions in many ways including 
impact parameter (R p ) distributions [8] and zenith angle distributions [11]. 
For this study, we presumed an average atmospheric clarity [12]. 

In order to calculate the autocorrelation function for this subset of data, 
we must first parameterize the HiRes-I monocular angular resolution. For a 
monocular air fluorescence detector, angular resolution consists of two com- 
ponents, the plane of reconstruction, that is the plane in which the shower is 
observed, and the angle if) within the plane of reconstruction (see figure 1). 
We can determine the plane of reconstruction very accurately. However, the 
value of if) is more difficult to determine accurately because it is dependent on 
the precise results of the profile-constrained fit [8,10]. 

The HiRes-I angular resolution is therefore described by an elliptical, two- 
dimensional Gaussian distribution with the two Gaussian parameters, and 
Cpiane, being defined by the two angular resolutions. For the range of estimated 
energies considered in this paper, = [4.9,6.1]° and <7 p i ane [0.4, 1.5]°. In fig- 
ure 2, the arrival directions of the HiRes-I events are plotted in equatorial 
coordinates along with their la error ellipses. 

In order to understand the systematic uncertainty in the angular resolution 
estimates, we consider a comparison of estimated arrival directions that suc- 



3 




Fig. 1. The geometry of reconstruction for a monocular air fluorescence detector 




Right Ascension (Hours) 



Fig. 2. The arrival directions of the HiRes-I monocular with reconstructed energies 
above 10 19 - 5 eV events and their la angular resolution 

cessfully reconstructed in both HiRes-I monocular mode and HiRes stereo 
mode. Because of the dearth of events with estimated energies above 10 19 - 5 eV 
that reconstructed satisfactorily in both stereo and mono mode, we consider 
all mono/stereo candidate events with estimated energies above 10 18 ' 5 eV. In 
stereo mode, the shower detector planes of the two detectors are intersected, 
thus the geometry is much more precisely known and the total angular reso- 
lution is of order 0.6°, a number that is largely correlated to o" p i ane and thus is 
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Fig. 3. Arrival direction error comparison between real data (mono vs. stereo) and 
simulated data for events with estimated energies above 10 18 ' 5 eV. The solid line 
histogram corresponds to the arrival direction error distribution of the monocular 
reconstructed Monte Carlo simulated data. The crosses correspond to the arrival 
directions error distribution observed for actual data by comparing the arrival di- 
rections estimated by the monocular and stereo reconstructions. The solid line in 
the ratio component corresponds to the fit y = ax + b where a = 0.000 ± 0.011 and 
b = 0.98 ±0.11. 

negligible when added in quadrature to the larger term, a^. This allows us to 
perform a comparison of the angular resolution estimated through simulations 
to the observed angular resolution values of actual data. In figure 3, we show 
the distribution of angular errors for real and simulated data. The uncertainty 
in the slope of the ratio (figure 3b) leads to an 7.5% uncertainty in the angular 
resolution. 





3 The Published AGASA Data 



The AGASA data with energies above 40 EeV has been published up to the 
year 2000 [2] and all but one of these events used for this calculation has a 
measured energy greater than 4 x 10 19 eV. The AGASA estimated angular 
errors [1] are shown in figure 4. The AGASA angular errors (figure 4) are fit 
to a two-component Gaussian distribution: 



n = N (E EeY ) 



O.33A#e-( A0 ) 2 / 2 ^ + O.Q7A6e~^ 2 / 2 ^ 



(1) 
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Fig. 4. The AGASA angular resolution as a function of estimated energy [1] 




Right Ascension (Hours) 

Fig. 5. The arrival directions of the published AGASA events with their 68% angular 
resolution 



where a x = 6.52° - 2.16° log 10 E EcV , a 2 = 3.25° - 1.22° log 10 E EcV , and N (E) 
is a numerically determined normalization constant. Figure 5 shows the arrival 
directions of the published AGASA events plotted in equatorial coordinates 
with their 68% angular resolution. 
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Fig. 6. An example of the autocorrelation function for a simulated data set that 
contains ~ 10 clusters in a total of 60 events — (a) the full autocorrelation func- 
tion for 9 = [0°,180°]; (b) the critical region of the the autocorrelation function: 
9 = [0°,10°]. 

4 The Autocorrelation Function 



We measure the degree of autocorrelation in both samples by means of an 
autocorrelation function. It is calculated as follows: 

(1) For each event, an arrival direction is sampled on a probabilistic basis 
from the error space defined by the angular resolution of the event. 

(2) The opening angle is measured between the arrival directions of a pair of 
events. 

(3) The cosine of the opening angle is then histogrammed. 

(4) The preceding steps are repeated until all possible pairs of the events are 
considered. 

(5) The preceding steps are repeated until the error space, in the arrival 
direction of each event, is thoroughly sampled. 

(6) The histogram is normalized and the resulting curve is the autocorrelation 
function. 

Figure 6a shows an example of the autocorrelation function for a highly clus- 
tered set of simulated data. The sharper the peak at cos # m i„ is, the more 
highly autocorrelated the data set is. There are many ways that one could 
quantify the degree of autocorrelation that a set possesses. The most obvi- 
ous way is to look at the value of the bin which contains cos6> min . However, 
this method has some pitfalls. First of all, the value of the last bin is depen- 
dent upon the chosen bin width. Also, the value of the last bin is not stable 
unless the angular resolution is sampled at a level that is computationally 
unfeasible. Finally, the value of the last bin over a large number of similarly 
autocorrelated sets does not produce a Gaussian distribution (see figure 7a), 
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Fig. 7. Distributions of normalized bin densities of cos0 m i n and <cos #>[o°,io°] val- 
ues for a large number of simulated sets with the same level of clustering as in fig- 
ure 6 — (a) Distribution of observed normalized bin densities of cos m i n , note that it 
is not Gaussian (x 2 /dof = 5.44); (b) : <cos #>[o°,io°] distribution (x 2 /dof = 1.09). 

thus complicating the interpretation of the results of an analysis employing 
cos# min as an observable. 

A more well-behaved measure of the autocorrelation of a specific set of data is 
the value of <cos 9> for 9 < 10°. This value is also a measure of the sharpness 
of the autocorrelation peak at cos 9 = 1. However, this method of quantifica- 
tion does not depend on bin width and it does produce Gaussian distributions 
when it is applied to large numbers of sets with similar degrees of autocorrela- 
tion as is demonstrated in figure 7b. An additional advantage to this method 
is that by considering the continuous autocorrelation function over a specified 
interval, both the peak at the smallest values of 9 and the corresponding sta- 
tistical deficit in the autocorrelation function at slightly higher values of 9 are 
taken into account. Thus we simultaneously measure both the positive and 
negative aspects of the autocorrelation signal. The interval of [0°, 10°] was 
chosen because in simulations it was found to optimize the autocorrelation 
signal for clusters resulting from point sources spread isotropically across the 
sky. 

Using the description of the HiRes-I monocular angular resolution above, we 
then calculate the autocorrelation function via the method described above. 
In figure 8, we show the result of this calculation. For this sample, we obtain 
<cos 0> [O o,io°]= 0.99234. 

We also calculate the autocorrelation function for the published AGASA 
events. We show the result in figure 9. For this sample, we obtain <cos #>[o°,io°] = 
0.99352. 
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(a) 

Fig. 8. The autocorrelation for the HiRes-I events above 10 19 ' 5 eV — (a) the full au- 
tocorrelation function for = [0°, 180°]; (b) the critical region of the autocorrelation 
function: <cos0>[ O o 1O o] = 0.99234. 




Fig. 9. The autocorrelation for the published AGASA events — (a) the full autocor- 
relation function for 9 = [0°,180°]; (b) the critical region of the autocorrelation 
function: <cos #>[o°,io°] = 0.99352. 

5 Quantifying the Relative Sensitivity of HiRes-I and AGASA to 
Autocorrelation 



In order to quantify the relative sensitivity of the AGASA and HiRes-I data 
sets, we must first understand the exposures of both detectors. For HiRes-I, 
we assemble a library of approximately 8 x 10 4 simulated events with energies 
above 10 19 ' 5 eV. We then pair each event with times during which the detector 
was operating. A mirror-by-mirror correction is applied where simulated events 
are rejected if the mirror(s) that would have observed the event in question was 
not operating at the time that event would have occurred. Once 10 7 pairings 
of simulated events and times are assembled, a surface plot is created of the 
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Fig. 10. Hires-1 estimated relative exposure, pn(6,a), for events above 10 19 5 eV in 
equatorial coordinates (right ascension right to left). The lightest region corresponds 
to a normalized event density of 2.5. The observable sky extends from 5 = —30° to 
5 = 90°. 

event density on a bin by bin basis. The value of each bin is then normalized 
so that the mean value of all the bins in the observable sky 8 = [—30°, 90°] is I. 
The resulting surface plot is shown in a Hammer-Aitoff projection in figure 10. 
We have previously shown that this method produced zenith angle, azimuthal 
angle, and sidereal time distributions that were consistent with that observed 
in the actual data [11]. The highest exposure areas have a normalized relative 
exposure: pu(5,a) =~ 2.5. 

For the AGASA detector, we refer to the distribution of event declinations 
presented in Uchiori et al. [13]. By following the lead of Evans et al. [14], we 
fit a normalized polynomial to this distribution: 

N(8) = 0.323616 + 0.03615155 - 5.04019 x HT 4 <5 2 + 

5.539141 x 10~ 7 5 3 ; (2) 

where N(5) holds for 5 = [—8°, 87.5°]the maximum value of N(5) is 1. We also 
know that: 

87.5° 87.5° 

A Q J N(5)d5= J p A (5) cos 6 dS, (3) 

-8° -8° 

where A Q is a numerically determined normalization constant. We then derive: 
p A (5) = A N(5) sec 5; A = 1.0251. (4) 
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Fig. 11. AGASA estimated relative exposure, pa(S), for events above 10 195 eV in 
equatorial coordinates (right ascension right to left). The lightest region corresponds 
to a normalized event density of ~ 1.6. The observable sky extends from 5 = —8° 
to 5 = 87.5°. 
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Fig. 12. Distribution of <cos 0>[o°,io°] values for simulated isotropic data sets — (a) 
HiRes-I; (b) AGASA. In each figure, the vertical line represents the the value of 
<cos 0>[o°,io°] f° r the- observed data. 

The value of each bin is once again normalized so that the mean value of 
all the bins in the observable sky 5 = [—8°, 87.5°] is 1. The resulting surface 
plot is shown in a Hammer-Aitoff projection of a equatorial coordinates in 
figure 11. The highest exposure areas have pa{ol) =~ 1.6. In figure 12, we 
show the distribution of <cos #>[o°,io°] values for isotropic data sets with each 
of the two different exposure models (HiRes-I and AGASA). The AGASA data 
set manifests ~ 10~ 3 chance probability above background. For the AGASA 
data, we also calculated the autocorrelation function without consideration 
to angular resolution and employed the more conventional 9 m i n observable. 
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After varying the bin width for 9 min and accounting for the trials factor, 
we independently concluded that the chance probability is ~ 10~ 3 for the 
optimal bin width, 9 min = [0°,2.5°]. We thus conclude that factoring angular 
resolution into our analysis and employing <cos 9> [o°,io°] as an observable in 
no way diminishes the sensitivity to autocorrelation in the reported AGASA 
data. 

There are a few important differences between the exposure of the HiRes- 
I and AGASA detectors. First of all, the exposure of the HiRes-I detector 
is more asymmetric than the exposure of the AGASA detector. This is not 
only due to seasonal variations in the HiRes detector, but also due to its 
ability to constantly observe the region around 5 = 90° due to a higher zenith 
angle acceptance. This higher zenith angle acceptance also allows the HiRes 
detector to observe a greater region of the southern hemisphere. In general, 
while AGASA reports observations for 56.9% of the total sky, the HiRes-I 
detector reports observations for 75% of the total sky. 

To simulate clustering we use the following prescription: 

(1) An event is chosen based upon the distribution in a and 5 that is dictated 
by p. In the case of HiRes-I, this is simply done by selecting a simulated 
event from our library and then assigning it a time that is a known good- 
weather ontime for the mirror(s) that observed that event. In the case 
of the AGASA detector, this is done by selecting a random value for S 
that conforms to the distribution in equation (4) and then assigning it 
a random value in a between Oh and 24h and sampling a value for the 
energy from the energies of the reported events. 

(2) This event does not represent the source location itself, but is assumed 
to have arrived from the source location with some error. We construct 
a "true" source location by sampling the error space of this event. 

(3) For each additional event assigned to that source, a simulated event is 
selected with a "true" arrival direction that is the same as that of the 
initial event. 

To study the relative sensitivity of AGASA and HiRes-I, we measure the value 
of <cos #> [00,10°] for multiple simulated sets with a variable number of doublets 
inserted. We then construct an interpolation of the mean value and standard 
deviation of <cos #>[o°,io°] from a given number of observed doublets for each 
experiment. This will allow us to state the number of doublets required for 
each experiment in order for the 90% confidence limit of <cos 9> [o°,io°] to be 
above the background value of 0.99250. Figure 13 shows the result of these 
simulations. In general, for the HiRes-I data set, the 90% confidence lower 
limit corresponds to the mean expected background signal with the inclu- 
sion of 6.25 doublets. For AGASA,the 90% confidence lower limit corresponds 
to the mean expected background signal with the inclusion of 5.5 doublets. 
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Fig. 13. Relative sensitivity of HiRes-I and AGASA to doublets — (a) Simulations 
with the HiRes-I detector and 52 events; (b) 90% confidence above background: 6.25 
doublets, 95% confidence above background: 8.25 doublets, (c) simulations with the 
AGASA detector and 59 events; (d) 90% confidence above background: 5.5 doublets, 
95% confidence above background: 7.0 doublets. In each figure, the horizontal line 
indicates the expected value of <cos #>[o°,io°] f° r an isotropic background 



This demonstrates that while AGASA has a slightly better ability to perceive 
autocorrelation, the sensitivity of the two experiments is comparable. 

We now apply the actual Hires-I <cos #>[o°,io°] to the sensitivity curve shown 
in figure 13. In figure 14 we can see the result of these simulations. The ob- 
served HiRes-I signal corresponds to the 90% confidence upper limit with the 
inclusion of only 3.5 doublets beyond random background coincidence. 

If we repeat this analysis with first, a 7.5% reduction in the estimated angular 
resolution values and second, a 7.5% increase in the estimated angular resolu- 
tion values, we obtain a range for the 90% confidence upper limit of [2.75, 4.0] 
doublets and a range for the 95% confidence upper limit of [4.5, 5.5] doublets. 



13 




0.991 




(a) 



5 10 

Number of Doublets 



0.991 



(b) 



5 10 

Number of Doublets 



Fig. 14. Sensitivity of the HiRes-I monocular observations to doublets — (a) Simula- 
tions with the HiRes-I detector and 52 events; (b) 90% confidence above observed 
signal: 3.5 doublets, 95% confidence above observed signal: 5 doublets. In each plot, 
the horizontal line represents the value of <cos #>[o°,io°] f° r the- observed HiRes-I 
data 

A final area of concern is the systematic uncertainty in the determination of 
atmospheric clarity. Because hourly atmospheric observations are not avail- 
able for the entire HiRes-I monocular data set, we have relied upon the use of 
an average atmospheric profile for the reconstruction of our data [12]. While 
different atmospheric conditions have negligible impact on the determination 
of the arrival direction for events with measured energies this high, differing 
conditions can have an impact on energy estimation and thus the number of 
events that are included in our data set. Over the la error space for our esti- 
mation of atmospheric conditions, the total number of events in our data set 
fluctuates on the interval [41,65]. The value of the observable, <cos#>ro°,io°], 
has a fluctuation on the interval [0.99226, 0.99249] owing to addition and sub- 
traction of events from the data set. Note that in neither case does the value 
of <cos 9> [o°,io°] exceed the mean value (0.99250) expected for a background 
set. 



6 Conclusion 



We conclude that the HiRes-I monocular detector sees no evidence of cluster- 
ing in its highest energy events. Furthermore, the HiRes-I monocular data has 
an intrinsic sensitivity to global autocorrelation such that we can claim at the 
90% confidence level that there can be no more than 3.5 doublets above that 
which would be expected by background coincidence in the HiRes-I monocu- 
lar data set above 10 19 ' 5 eV. From this result, we can then derive, with a 90% 
confidence level, that no more than 13% of the observed HiRes-I events could 
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be sharing common arrival directions. This data set is comparable to the sen- 
sitivity of the reported AGASA data set if one assumes that there is indeed a 
30% energy scale difference between the two experiments. It should be empha- 
sized that this conclusion pertains only to point sources of the sort claimed by 
the AGASA collaboration. Furthermore, because a measure of autocorrelation 
makes no assumption of the underlying astrophysical mechanism that results 
in clustering phenomena, we cannot claim that the HiRes monocular analysis 
and the AGASA analysis are inconsistent beyond a specified confidence level. 
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